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We propose an experimental setup using 3D Dirac semimetals to access a novel phenomenon 
induced by the chiral anomaly. We show that the combination of a magnetic field and a circularly 
polarized laser induces a finite charge density with an accompanying axial current. This is because 
the circularly polarized laser breaks time-reversal symmetry and the Dirac point splits into two 
Weyl points, which results in an axial-vector field. We elucidate the appearance of the axial-vector 
field with the help of the Floquet theory by deriving an effective Hamiltonian for high-frequency 
electric fields. This anomalous charge density, i.e. the chiral pumping effect, is a phenomenon 
reminiscent of the chiral magnetic effect with a chiral chemical potential. We explicitly compute 
the pumped density and the axial-current expectation value. We also take account of coupling to 
the chiral magnetic effect to calculate a balanced distribution of charge and chirality in a material 
that behaves as a chiral battery. 


I. INTRODUCTION 

Quantum anomaly^’^ is one of the most fundamental 
concepts in quantum field theory and anomaly-related 
phenomena appear in various physical systems from ele¬ 
mentary particle reactions at high energy to tabletop ex¬ 
periments in condensed matter physics. Thus, it would 
offer tremendous opportunities in theory and experiment 
to seek for novel manifestation of quantum anomaly. 

Recently, the chiral magnetic effect (CME) has been 
attracting broad interest, which was first triggered in 
the research field of the ultra-relativistic heavy-ion col¬ 
lision^. Later on, the CME theory has been formulated 
in terms of the chiral (axial) chemical potential /ta"* and 
in this way the application opportunities of the CME 
have opened to a wider range of the research fields. Al¬ 
though it is still challenging to confirm any smoking- 
gun experiment through the charge asymmetry fluctu¬ 
ations in the heavy-ion collision®, materials in condensed 
matter systems would provide us with more controllable 
environments. Recently, the material realization of 3D 
Dirac fermions®^^ and Weyl fermionshas been es¬ 
tablished. In fact, it is claimed that an evidence for the 
CME has been observed in ZrTes with parallel electric 
and magnetic fields, for which the magnetoconductance 
has quadratic dependence on the magnetic field^'^. 

In the present paper, we propose another manifestation 
of quantum anomaly which we call the “chiral pumping 
effect” (CPE). This effect can be understood as a cousin 
of the CME originating from a variant of anomaly rela¬ 
tion. First, let us recall that the CME relation^ is given 
by a compact formula: 

f = oc (CME) , (I) 

27r^ 

where represents the Levi-Civita symbol. This rela¬ 


tion states that a charge current is induced in parallel 
to the applied magnetic field provided that 

the chiral chemical potential (= /ta) is nonzero. We 
can anticipate the CPE by exchanging the 3 and 0 in¬ 
dices in the above relation, which leads to the following 
relation for a charge density: 

f oc B (CPE) . (2) 

Of course, the net charge must be conserved, so we must 
understand this relation in an open system setup. This 
means that the charge is pumped into the system from 
the surrounding reservoir, and the change is proportional 
to the inner product of the axial-vector field A 5 and 
the magnetic field. In the language of Weyl semimetals, 
A 5 corresponds to the displacement of two Weyl points. 
Since the CPE originates from quantum anomaly simi¬ 
larly to the CME, we expect that the coefficient in the 
relation ( 2 ) should be renormalization free. 

Here, let us mention on a related approach with an 
inhomogeneous axion term 9{x, t)E ■ B' The Chern- 
Simons-Maxwell equations suggest that a charge density, 
jo oc (V0) • B, is induced^^. On the algebraic level this is 
almost equivalent to Eq. (2) once V0 is identified with A 5 
(which is possible only in the massless case). The crucial 
difference in physics is that 9{x,t) should usually belong 
to the property of the material, whereas our A 5 in Eq. (2) 
is an external field that is experimentally adjustable as 
we see later. It is actually quite non-trivial how to im¬ 
plement such axial-vector field using electromagnetic de¬ 
vices. In the case of experimental confirmation of the 
CME, in fact, it is technically difficult to control ^a and 
so the formula ( 1 ) cannot work as it appears for the ex¬ 
perimental signature. It should be noted that 9(t) = 
in the Chern-Simons-Maxwell theory would immediately 
lead to Eq. (I) as pointed out originally in the very first 
CME paperand later in the context of Weyl semimetals 
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Circularly Polarized Laser 



FIG. 1. Experimental setup with the magnetic field and the 
circularly polarized laser onto a 3D Dirac semimetal. 


also’^®. In contrast, our case of the CPE has an advan¬ 
tage that we can easily manipulate A^. Moreover, the 
balanced conhguration of charge and axial-charge (i.e. 
chirality) turns out to be a system of capacitor of chiral¬ 
ity which should be useful for more direct CME studies. 

The aim of this work is to propose a tractable experi¬ 
mental setup to manifest the CPE in 3D Dirac systems. 
A key step to realize the axial-vector field experi¬ 
mentally is, as discussed below, that we utilize a rotating 
electric field, i.e., circularly polarized laser rotating in a 
two dimensional plane (see Fig. 1 for a schematic illus¬ 
tration). We also refer to a related idea with circular 
polarizations in 3D Dirac semimetals^® and more general 
photo-induced effects^'’. Using a simple fermionic de¬ 
scription, we will show that the Dirac point splits into two 
Weyl points. With an additional magnetic field^^, a finite 
density arises from the lowest Landau level (LLL) of one 
chirality, which manifests a concrete picture of the CPE 
in (l-l-l)-dimensionally reduced theory of fermions^^. 

This paper is organized as follows. In Sec. II we discuss 
the Floquet effective Hamiltonian to conhrm an axial- 
vector field. In Sec. Ill we consider a combination with 
a magnetic field and perform explicit calculations for the 
charge density and the axial current. Inhomogeneous 
electric charge and chirality should be balanced with each 
other. We solve these coupled equations of the CPE and 
the CME to obtain a balanced distribution of the elec¬ 
tric charge and the chirality in Sec. IV. Finally, Sec. V is 
devoted to our discussions and conclusions. 


II. FLOQUET EFFECTIVE HAMILTONIAN 
AND AXIAL-VECTOR FIELD 

We explain how to realize the axial-vector held in a 
3D gapless Dirac system by applying a circularly polar¬ 
ized laser. We note that concrete calculations below are 
known ones, but a clear recognition of the axial-vector 
held has not been established. When continuous laser 
helds are imposed externally, the Hamiltonian H{t) be¬ 
comes periodically time-dependent, i.e., H{t + T) = H{t) 
where T = 27:10. is the periodicity. Quantum states in 
time periodic driving are described by the Floquet the- 
ory2-'^>24^ that is, a temporal version of the Bloch theorem. 


The essence of the Floquet theory is a mapping between 
the time-dependent Schrodinger equation and a static 
eigenvalue problem. The eigenvalue is called the Floquet 
pseudo-energy and plays a role similar to the energy in 
a static system. Applications of the Floquet theory to 
periodically driven systems with topology changing has 
been a recent hot topic^®^^® and experiments have also 
been done^®^®^. 

To make this paper as self-contained as possible, in this 
section, we derive the effective Hamiltonian in an explicit 
way, though the final result is not very new but already 
known. Let us consider a Hamiltonian, iLtot = Hq Hmt, 
with 


Ho = 7°7 • p + 7 °m , Hint = -e7°7 A , (3) 


that describes the one-particle Dirac system coupled to 
an external gauge field and are the Dirac matrices 
satisfying { 7 '', 7 '^} = 2^^^. In an electric field with cir¬ 
cular polarization in the x-y plane, we can write the time- 
dependent vector potentials down as 

^ cos(Qt) , Ay = ^ sin(Dt) A^ = 0 , (4) 


where O is the frequency. We can conveniently de¬ 
compose the interaction part of the Hamiltonian into 
two pieces as Hint = -I- where H± = 

— with 7 ^ = 5 ( 7 “ ± ^ 7 ^). Now we assume 

that the the period T = 27r/D of the circular polarization 
is small enough as compared to the typical observation 
timescale. We can then expand the theory in terms of 
Lo/O (with ui being a frequency corresponding to some 
excitation energy). Taking the average over T we can 
readily find the following effective Hamiltonian^®’®®^®®: 


In 




c^Ho + -[H_,H+] , 


( 5 ) 


to the first order in the expansion. We can also find 
the same form from the Floquet Hamiltonian using the 
Van Vleck perturbation theory®®. Interestingly we can 
express the induced term as 


Hino = 1^\H_,H+] = , ( 6 ) 


where we defined /3 = (eA)^/D®. This means that the 
circular polarized electric field would induce an axial- 
vector background field A 5 = Pz perpendicular to the 
polarization plane. Essentially the same expression was 
obtained in the context of “Floquet Weyl semimetal” and 
the corresponding Floquet bands were hgured out®®. The 
physics is basically the same as this preceding work®®, but 
we use a different language here and, for completeness, 
we shall see the energy dispersion relations in the rest of 
this section. 

The effect of finite /3 is easily understandable from the 
energy dispersion relations. We can immediately diago¬ 
nalize iLeff and the four pseudo-energies read: 


e±{p) = yJpl+pl + {^/pl+m‘^±py 


( 7 ) 
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FIG. 2. Pseudo-energies ±£-(p) as a function of = 
\/p| + Py fl'iid for /3 > m [in the fignre {m, fi) = (1,5) 
was chosen]. 


and —s±{p). We display low-lying ±S-{p) in Fig. 2. 
which shows that the Dirac point splits into two Weyl 
points with a displacement given by 

Ap = y/— rn? . ( 8 ) 

In fact, P is nothing but a momentum shift along the 
z-axis that is positive for the right chirality state (i.e., 
75 '0 R = +V'r) and negative for the left chirality state 
(i.e., ysV'L = —'0 l)- We point out that time-resolved 
ARPES should be able to see this splitting of Weyl points 
in a similar manner as the gap opening"’ of the 2D Dirac 
point already observed experimentally^ b Interestingly, 
as long as /? > m, the pseudo-energy always has two 
Weyl points (if they are inside of the Brillouin zone) 
even for m > 0. Therefore, we do not have to require 
strict massless-ness to realize gapless dispersions, which 
should be a quite useful feature for practical applications 
including the Schwinger or Landau-Zener effect. In what 
follows below, we limit ourselves to the m = 0 case just 
for technical simplicity. 

The generalization from the one-particle Hamiltonian 
to the many-body field theory is straightforward. It is 
then more convenient to work with the Lagrangian den¬ 
sity corresponding to H^g, that we can express as 

Ceff = - m)ip + /3'07^75V' + /rA'07°75'0 • (9) 

Here we include pA for completeness, which is a necessary 
ingredient for the CME. It is clear from this Hamiltonian 
that we should identify /3 as a parameter representing 
what is called the chiral . We should emphasize 

a crucial difference from the idea of the chiral shift that 
is not directly controllable but secondarily generated by 
finite-density effects. In our present problem, however, 
P is an external parameter that we can control with the 
amplitude and/or the frequency of the circular polarized 
electric field. What we will see is that, conversely to 
discussions on the chiral shift^^’^®, a finite density is gen¬ 
erated by this externally given P ^ 0. 



FIG. 3. Energy dispersion relations of the LLL (i.e., = 

= 0) with (a) chiral chemical potential pA relevant for 
the CME and (b) chiral shift /3 relevant for the CPE. Filled 
(blank) dots represent states that are newly occupied (unoc¬ 
cupied) . 


III. RESPONSE TO THE MAGNETIC FIELD 

It is the most essential point that we can regard P 
as the z component of an axial-vector field; /3 ~ A|. 
Then, if we further impose an external magnetic field 
on this system, we should expect the following 
anomaly relation; oc PB, immediately from the CPE 
or Eq. (2). We will confirm this expectation with explicit 
calculations, but before going into details, let us consider 
an intuitive interpretation to understand Eq. (2). 

For the purpose of comprehensible illustration it would 
be useful to sketch the dispersion relations in the same 
way as in a CME literature'*. Figure 3 (a) shows the 
dispersion relations of the lowest Landau levels (LLLs) 
relevant for the CME with P = 0 and pA ^ 0. In this 
case with pA > 0 the energies of the right-handed (i?) 
particles are decreased, while those of the left-handed 
(L) particles are increased. Note that only one spin state 
is chosen out for the LLL depending on the sign of eB. 
Therefore, a positive pA favors more R than L. This 
explains how a finite chiral density is accumulated, while 
a net density remains vanishing. Also we see that the 
CME current flows from the LLL with pA and gives an 
intuitive picture for the anomaly relation (I). In contrast, 
as seen in Fig. 3 (b) for p > 0, the energy dispersion of R 
(and L) is shifted positively (and negatively, respectively) 
along the p^ axis. Thus, assuming that particles can flow 
in through the p^-integration edges [as indicated by the 
dashed lines in Fig. 3 (b)], both L and R LLL states 
increase their occupation. This is the mechanism of how 
a finite density develops from a combination of P and 
B (that causes the dimensional reduction) through the 
anomaly relation (2). The argument at the same time 
tells us that a R excess in pz < 0 and an L excess in 
Pz > 0 contribute to a negative flow of the net chirality, 
which is an analogue of the CME current. Summarizing 
the above discussions we should expect non-zero (j°) and 
01 )- 

Now we shall make sure of our qualitative argument 
by performing explicit calculations and identify the co¬ 
efficient of Eq. (2) following the calculations done previ¬ 
ously'**"’®®. We set Pa = 0 in the following and concen¬ 
trate on the CPE only. The fermion propagator with P 
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in the presence of eB > 0 is 

OO 

G{x ^0)^ Go-(a;)P_ + ^[G+(a:)P++G-(x)P_] (10) 

n—1 


near the coincidence limit, where P± = |(1 ± i 7 ^ 7 ^) is 
the spin projection operator and the first term, GqP-, 
is the LLL contribution. The (Feynman) propagator for 
each Landau mode is 


G^(x^O) 


2n 


duj dp 


Kn{uj,±{p + I3))P^ 


+ Kn{uj,T{p- I3))P5 


( 11 ) 


Here, P^ = 1(1 ± 75 ) is the chirality projection operator 
and we defined 


Kn{uj,p) 


LO +p 

Up — 2eBn — p'^ + ie 


( 12 ) 


Using this form of the propagator we can express the 
current expectation value as 


if) 


-e lim tr[ 7 °G(t)] . 

t-s-0+ 


(13) 


We can simplify the calculation by noting that 
tr(Pg^P±) = 3 tr[(l ± 75)(1 ± = 1 where two ± 

are independent. Because lu in the numerator should be 
vanishing after the w integration, it is straightforward to 
check that tr( 7 °G+P+) + tr( 7 °G“P_) = 0. Therefore 
only the LLL contribution survives, which yields: 


e^B f f PP P P — P \ 

2{2tt) J 27r\ \p +\p-P\) 


(14) 


after the uj integration. The above is a finite integral 
only from —fd < Pz < P- Actually, this simple expression 
is precisely a concrete realization of our qualitative argu¬ 
ment based on Fig. 3 (b). Finally, after the p integration, 
we recover the same result as previously obtained one’^^: 


(J°) = 


e^PB 

~2^ ’ 


(15) 


which is quite reminiscent of the formula for the CME 
current, namely, (j^) = /{2tt‘^). This compact for¬ 

mula itself is a known one, but the fact that the above 
(j^) can be easily realized with a circularly polarized laser 
is our main claim in this work. 

Remembering our qualitative argument, we can in turn 
expect an axial current (j^) which is similar to what is 
called the chiral separation effect; (j^) oc pB^^. To see 
this, let us perform an explicit calculation using (jjP) = 
—etr[ 7 ^ 75 G( 0 +)]. In this case we use tr( 7 ^ 757 °P^P±) = 
itr[7"757°(l±75)(l±*7"7^)] = ±ttr(7^757V7^) = 
±1, where the last ± refers to ± of P±. Then, the LLL 
contribution is completely identical to that for (j°) ex¬ 
cept for the overall sign, that means, (iA)(LLL) = 

We can immediately explain this proportionality for the 


LLL from the property of the two-dimensional Dirac ma¬ 
trices; 7 '^ 75 ^^ = 7 ° where 75 ^^ is the two-dimensional 
counterpart defined by 7 ^ 7 ^. In this case, however, some 
complication appears from non-zero Landau levels as is 
the case for the chiral shift'^^’^®. We can write those con¬ 
tributions down explicitly as 


Ua) (n> 0 ) 


e^B 

2'k 



X - 


P + P 


P-P 


yj2eBn -I- {p+PY y/2eBn + {p—PY 


(16) 


This is a subtle expression whose precise value depends 
on how to organize the infinity in the Landau sum and 
the momentum integration. For example, if we take 
the Landau sum up to N and the momentum inte¬ 
gration within \p\ < A in such a way that 2eBN 
A^, then the p integration is easy to perform and we 
find {jA)(n>o) = ~2Af(j°). For more general situation 
N should be replaced with {A/eB){y/2eBN -|- A^ — A), 
which approaches yj2NhP/eB in the opposite limit of 
A^ ^ 2eBN. In reality these cutoffs should be hxed 
by the microscopic properties of the material, especially 
the Debye mass and the Brillouin zone structures. The 
conclusion is that the system comes to have a finite axial- 
current as given by 


/■z\ e'^PB 

(■?a) = -« , (17) 

the coefficient a of which is not anomaly protected unlike 
Eq. (15) and it is a material-dependent problem to fix a. 


IV. STATIC CHARGE DISTRIBUTION 

So far, we assumed an infinitely large system. Impos¬ 
ing surface boundary conditions with a finite extent, we 
should take account of the polarization and screening ef¬ 
fects to find a balanced distribution of the charge and 
the chirality density. Let us consider a finite-size mate¬ 
rial whose thickness in the z direction is d; we take the 
z coordinate so that the material is placed in a range 
—d/2 < z < -l-d/2. Then, because of the surface effects, 
we should introduce z dependent p{z) and PAiz), the 
determination of which is the goal of this section. 

In the presence of pa{z), the CME current should fol¬ 
low from Eq. (1). In equilibrium, however, there should 
be no current and the CME current should be canceled 
by the polarization effect, i.e. 

/total = + ^/ dzf{z) , (18) 

where the second term in the right-hand side appears 
from Ohm’s law with the electric conductivity ct > 0. In 
fact, J dz j°(z) is nothing but an electric field associated 
with the polarization or the charge distribution j®(z). 
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Qualitatively, this last term represents a flow to flatten 
the charge distribution. For the balance equation for 
we postulate the same structure as Eq. (18) with the 
replacement of and ^ (i-e., the chiral 

separation effect'^'*): 

Ja .total = . (19) 

where we added the last term to consider diffusion pro¬ 
cesses of chirality introducing a diffusion constant A > 0. 
Intuitively, this last term represents movement of chiral¬ 
ity to decrease the chirality gradient. The densities, 
and j^, are given as 


e'^B ^ e , , o 


•’» = to^***'*)’• (*> 


where a CPE contribution to is added to the stan¬ 
dard density-chemical potential relation. We plug these 
expressions to the conservation laws; = 0 and 

= 9 (note that there is no El • S with our elec¬ 
tromagnetic configuration). Since we are interested in 
the static profile only, we can drop the time-derivatives 
to reach finally: 


|^9,/rA(^)+a|^^+^M^)3 = 0, (21) 

^9^/i(z) - = 0 . (22) 


We should solve these differential equations under the 
charge constraint; 2 ^{z) = (e^il/27r^)/3c?, which 

leads to fiA{df2) — fiA{—d/2) = —/dad. If the thickness 
d of the material is sufficiently small, we can find an 
approximate solution by expanding the solutions in terms 
of z, that is. 


Kz) cs ^ , fiA{z) -Pcrz . (23) 

eB 

Here, we further used a condition to impose zero net chi¬ 
rality. It should noted that these are leading-order results 
and, to satisfy Eq. (21) strictly for example, we should 
consider a term oc zd in ^aIz) which we neglected. The 
charge density is thus more screened near the surface 
at larger \z\ unlike an ordinary conductor, while there 
are more R (and L) near the bottom at negatively large 
z (and near the top at positively large z, respectively). 
This is very interesting result; the system behaves as a 
capacitor of chirality or the chiral battery in which a non¬ 
zero slope in the chirality distribution is sustained by the 
CPE as well as a finite net density. Such a realization 
of the chiral battery should be useful for more quantita¬ 
tive investigations of the CME and related phenomena, 
e.g., the chiral plasma instability^'* for instance. We note 
that the idea of the chiral battery can be traced back 
to the original work and extensively discussed as chiral 
electronics ' . 


V. DISCUSSIONS AND SUMMARY 

The anomalously induced charge density (j**) given by 
Eq. (2) or (15) is an experimentally detectable quantity. 
One possible way to observe this is using transient grating 
techniques'*^. The grating of the circular polarized field 
would lead to inhomogeneous charge distribution, and 
the dynamics afterward can be measured. 

It is also an interesting possibility to make use of 
Eq. (17) instead of Eq. (15) since a could take a large 
number (involving some cutoffs at least in a naive es¬ 
timate). In the two-component spinor representation, 
Ja — is nothing but the spin expec¬ 

tation value. Hence, for systems with gapless Dirac par¬ 
ticles, Eq. (17) indicates not only the dynamical flow of 
chirality but also the spin polarization. Spin polariza¬ 
tion can be probed by pump-probe magneto-optical Kerr 
effect. In addition, when topological current or density 
exists in general, there can be some photo-emission pro¬ 
cesses via anomaly as discussed in Refs. '*'* ’'*“*. It deserves 
further investigations to quantify emitted photon spectra 
associated with the CME and the CPE. We will leave 
them for interesting future problems. 

In summary, we formulated the CPE that generates a 
finite density from the combination of the circular polar¬ 
ized electric field and the magnetic field. The interesting 
point in this setup is that the externally imposed elec¬ 
tric and magnetic fields are always pointing orthogonal to 
each other and there is no inner product oi E ■ B which is 
usually the source of topological charge and parity break¬ 
ing. Instead, the circular polarization breaks the parity 
symmetry, and so it alone can take care of the role of 
E ■ B. Besides, thanks to rotating electric field, we do 
not need to subtract a huge background due to ordinary 
electric current that obeys Ohm’s law. 

We gave intuitive arguments about the generation of 
the density and the axial current and confirmed our ex¬ 
pectations by explicit calculations recovering previously 
known expressions. We saw that the density (j**) is ex¬ 
pressed in a compact form similar to the CME and its 
coefficient is anomaly protected. In contrast to this, (j^) 
has complicated contributions that need some ultravio¬ 
let regularization or physical cutoff. We would here em¬ 
phasize the following point: Each building block for our 
conclusion was known; the appearance of a term 

from the Floquet Weyl semimetal and the topologically 
induced density from the chiral shift were all “re-derived” 
in this work, but our novelty is found in a combination of 
them. The most important is that our controlling param¬ 
eter /3 is given externally unlike the axion term associated 
with intrinsic properties of the material. 

Although it is difficult to determine the coefficient of 
(ja), it is still proportional to /3 or a combination of 
{eE)‘^/il^ (as long as H is large enough to justify the 
leading-order expansion as we did in this work). This 
characteristic dependence of should give a consis¬ 
tency check for anomalously induced (_)**) and {Ja)- 

Fortunately, such ambiguity in (j^) does not enter the 
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determination of the static charge and chirality distribu¬ 
tion in the material. Our study implies that the induced 
(jO) gg^g slightly smaller near the material surfaces. The 
most interesting is that the chirality has a non-trivial 
distribution also and the chirality separation is sustained 
by the CPE, which realizes a system of the chiral bat¬ 
tery. The experimental confirmation of the CPE itself 
is quite challenging, and furthermore, the CPE opens a 
new possibility for more direct CME studies using this 
chiral battery. 
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